Griffiths singularities in unbinding of strongly disordered polymers 
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Griffiths singularities occurring in the unbinding of strongly disordered heteropolymers are stud- 
ied. A model with two randomly distributed binding energies —1 and — v, is introduced and studied 
analytically by analyzing the Lee- Yang zeros of the partition sum. It is demonstrated that in the 
limit v — ► oo the model exhibits a Griffiths type singularity at a temperature To = 0(1) corre- 
sponding to melting of long homogeneous domains of the low binding energy. For finite »>1 the 
model is expected to exhibit an additional, unbinding, transition at a high temperature Tm = 0(v). 
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The unbinding transition of two polymers has at- 
tracted considerable attention in recent years. This phe- 
nomenon arises in a variety of physical contexts in which 
two linearly extended objects interact. Examples are 
the denaturation transition of double stranded DNA |l]] 
where the two strands of the molecule unbind; wetting in 
two dimensions 0] where the linear interface separating 
two bulk phases becomes unbound from the substrate; 
the depinning of flux lines from a columnar defect in 
type II superconductors Q and others. In these prob- 
lems an attractive interaction binds the two polymers, 
or the other linear objects, to each other at low tem- 
peratures. However, as the temperature is increased to 
a critical temperature a transition takes place at which 
they become unbound. 

In the case of homopolymers the unbinding transition 
is fairly well characterized. Using either a directed poly- 
mer approach Q or the more general Poland-Scheraga 
(PS) type models || the nature of the transition and the 
critical behavior have been studied. In the PS approach, 
introduced in the context of DNA denaturation, a mi- 
croscopic configuration of the two polymers is viewed as 
an alternating sequence of bound segments and denat- 
urated loops. A statistical weight is assigned to each 
bound segment and open loop. The weight of a bound 
segment of length I is taken as w with w = exp(— j3e), 
(3 = l/(ksT) being the inverse temperature, and e is the 
binding energy of two monomers. On the other hand a 
loop of length I is assumed to carry no energy but rather 
a degeneracy factor which takes the form O(i) ~ s l /l c . 
Here s is a geometrical non-universal constant while c is 
a universal constant. 

The nature of the unbinding transition was found to 
depend only on the parameter c Q. The transition is 
continuous for 1 < c < 2 and is first order for c > 2. No 
phase transition takes place for c < 1 and the polymers 
are bound at all temperatures. The value of c taken in 
the PS approach is universal and depends only on the 
dimension d of the system and on wether or not long 
range interactions, such as self-avoiding interactions, are 
taken into account [|| [| [?]. In the directed polymer- 
approach in d+ 1 dimensions c is related to the number of 



random walks which return to the origin for the first time. 
This is known to be given by || c = 2 — d/2 for d < 2, 
c = d/2 for d > 2. In d = 2 one has Q(l) ~ s l /(l)n 2 I). 

In many cases one is interested in the unbinding tran- 
sition of two heteroploymers which are composed of a 
sequence of two or more types of monomers. This is the 
case, for example, in DNA where each strand is composed 
of a sequence of four types of nucleotides. Heteropolymcr 
unbinding is usually modelled by introducing a quenched 
disorder in the contact energy between the two chains. 
This disorder, besides modifying the nature of the un- 
binding transition, may introduce Griffiths singularities 
in the free energy as is the case in, for example, random 
magnets ||. In this Letter we study these singularities 
in strongly disordered heteropolymer unbinding. 

For weak disorder the Harris criterion jl(J may be ap- 
plied in order to determine wether or not the critical 
behavior of the homopolymer unbinding is modified by 
the disorder. One finds that weak disorder is irrelevant 
as long as c < 3/2 [|ll] [l2| [l3| . This corresponds to 
1 < d < 3 for directed polymers. The marginal c = 3/2 
case, corresponding to d = 1 and d — 3 in directed poly- 
mers, has been studied in detail ||l2|, |l3|, 14}, 1^, |l6| and 



is found to be marginally relevant |L3|, |14 , |15 1 . 

Recently the case of strong disorder has been consid- 
ered by Tang and Chate Q . In this study the unbinding 
transition is analyzed using a real-space renormalization 
group scheme and numerical studies using a transfer ma- 
trix approach. It is suggested that in the strong disorder 
regime the unbinding transition becomes very smooth, of 
a Kosterlitz-Thouless (KT) type, with a weak free energy 
singularity of the form exp(— a/ \J\T — Tm\) for T < Tm- 
The free energy is differentiable to all orders at T = Tm- 
On the other hand, an analysis of the transition by Azbel 
p8[ suggests that the transition could be of very high but 
finite order (over 100). The order of the transition de- 
pends on non-universal quantities like the details of the 
distribution of the binding energies. 

In this Letter we introduce a simple model for the un- 
binding process of heteropolymers in the strong disorder 
regime. We demonstrate that the free energy exhibits a 
Griffiths type singularity @, [HI [2fl ^l| at a low tempera- 
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ture, Tq, resulting from the melting of long homogeneous 
domains of sites with low binding energy. The model then 
exhibits an unbinding transition at a higher temperature, 
T M - 

The model considered in this study is a PS model with 
two types of randomly distributed binding sites. The 
energy ej of site i is given by: 



— 1 , with probability p 
—v , with probability 1 — p 



(1) 



We first study the limit v — > oo, where no unbinding 
takes place. We demonstrate that a Griffiths singularity 
occurs at a temperature Tq = 0(1). We then argue that 
the transition persists at finite but large v, where it is 
followed by an unbinding transition at a higher tempera- 
ture Tm = 0(v). This model is similar in spirit to a class 
of models used in the past to study, for example, random 
field fsing models p2| . 

The simplifying feature of the limit v — > oo is that the 
chain is decomposed into independent non-interacting fi- 
nite segments. The free energy of the model is thus ob- 
tained by properly averaging the free energy of finite seg- 
ments, each of which is homogeneous. The free energy 
density of the chain may thus be expressed as 



OO 

N=0 



(2) 



Here t = (Tq — T)/Tq, where Tq is the melting temper- 
ature of a homopolymer with e — —1. The free energy 
density /jv(t) is that of a homogeneous chain of length 
N + 2 where the first and last sites of the two stands 
are taken to be bound. In the thermodynamic limit the 
singular part of the free energy is given by 



Jim f N (t) 

iv— »oo 




(3) 



where v is the critical exponent associated with the di- 
vergence of the correlation length £ ~ t~ v ' . The exponent 
is given by [§ 



l/(c- I) , for 1 < c < 2 



for c > 2 . 



(4) 



Clearly, at any finite N, /jv(i) is an analytic function 
of t and the singularity develops only in the limit N — > 
oo. The sum is a canonical case where Griffiths type 
singularities may arise. 

We proceed by first demonstrating that the average 
free energy (||) is singular at t = 0. This is done by 
showing that t = is an accumulation point of the poles 
of the free energy in the complex t plane. The average 



free energy is then proved to be differentiable to all orders 
at the transition. 

To show that the free energy (0) is singular at t = 0, 
we use the Lee- Yang theory and analyze Jm (t) in the 
complex w — exp(/3) plane for large N. The partition 
sum Zn(w) of a finite chain is a polynomial of order N 
or less in the complex w — wr + iwi plane. It has no 
zeros on the positive real w axis. Let pn{wr,wi) be the 
density of zeros in the complex w plane. The free energy 
density /n(w) may be expressed as 

f N (w) = -j^ J p(x, y) \n(w - z)dxdy (5) 

where z = x + iy. To calculate p(x, y) we note that it 
can be expressed as 23, p3| 



p{x,y) = 2ttV 2 



with 



lim 



ln|Z 



N 



(6) 



(7) 



One may thus view —<f> as an electrostatic potential in 
a two-dimensional space generated by a charge density 
p(x,y). We proceed by calculating near the transition 
where w — wq oc t. It is evident from Eq. (|3|) that the 
system exhibits two phases: an unbound phase (I) near 
the real negative t axis where = 1 and a bound phase 
(II) near the real positive t axis where = exp(Nt u ). 
In the complex t plane the potential in the two phases is 







, phase I 



r v cos(9u) , phase II , 



(8) 



where t = rexp(i8). The boundary between the two 
phases is found by requiring that </> is continuous. This 
yields the lines 9 — ±7r/ (2v). Hence phase II exists inside 
the wedge — tt/(2v) < 9 < +n/(2v) while phase I exists 
on the outside (see Fig. [I]). The discontinuity in the 
field E = \7(j) in the direction perpendicular to the phase 
boundary yields the charge density and thus the density 
of zeros of the partition sum on that line. The angular 
component of the electric field Eg — l/rdgcf) in phase II 
on the transition lines 9 = ±7r/(2f) is given by 



En 



(9) 



Thus, the density of the zeros of the partition sum along 
these lines in the N — > oo limit is 



p{r) 



2tt 



(10) 



It is evident that for c > 2, where v = 1, the density 
of zeros is constant. Thus the typical distance ro to the 
first zero is of order ro ~ 1/N. For I < c < 2, where 
v = l/(c— 1), the density of zeros vanishes as t — ► 0. The 
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FIG. 1: A schematic illustration of the lines (in gray) in the 
complex t plane on which the Lee- Yang zeros of the partition 
sum lie. The density of zeros on these lines satisfies p ~ r" _1 . 



FIG. 2: The Lee- Yang zeros of the partition sum ( |ll[ ) for a 
direct walk model corresponding to c = 3/2 are shown for 
N = 100 (inner points) and N = 380 (outer points). 



typical distance, r , is determined by r v ~ x dr ~ 1/JV. 
This yields ro ~ 1/N C . Note that this distance goes 
to zero in the thermodynamic limit as long as c > 1. 
On the other hand for c < 1, ro does not vanish in the 
thermodynamic limit and thus no singularity develops 
in the partition sum and no transition takes place, as 
expected. 

We have checked this picture by numerically evaluating 
the zeros of the partition sum of a simple model of a 
directed polymer in 1 + 1 dimension corresponding to c = 
3/2. In this model we consider a random walk which at 
each step moves either up (x — > x+1) or down {x — > x—1) 
with probability 1/2. We consider walks of length N such 
that xq = xn — for which x > at all intermediate 
steps. With each visit to the origin, x = 0, we associate 
a weight w. The partition sum of this model is 



N/2 

Z N = T w R ' 



R=l 



(11) 



where Qn.r — ( jy/2 ) ^ s ^ ne num b er of such walks 
of length TV which return to the origin R times |25| . It 
is straightforward to show that in the N — > oo limit, the 
model undergoes a transition at w = 2. The zeros of the 
partition sum ( |ll| ) have been calculated numerically for 
N = 100 and TV = 380 and are shown in Fig. |. It is 
readily seen that the zeros are located on a curve in the 
complex w plane. As N increases this curve develops a 
wedge near w = 2. In the large N limit the curve should 
approach the point w — 2 at an angle 9 = ±7r/4. 

We now turn to the average free energy density J 7 , Eq. 
|2|. It is clear that t = is an accumulation point of the 
poles of this free energy (or cquivalently of the zeros of the 
corresponding partition sum). Thus the series expansion 
of the free energy at t — has a zero radius of convergence 
and the free energy is singular. However, this singularity 



is expected to be rather weak since the weight of the 
poles located at a distance O^Y/N 1 /") from t = (which 
correspond to fisr(t) hi (|J)) have an exponentially small 
htp N . This is a typical case of a Griffiths singularity. 
To demonstrate that the average free energy J-(t) is 
differentiable to all order at t = we adopt the approach 
of |H] used for analyzing the Griffiths singularities in 
random ferromagnets. We start by considering the mth 
derivative of fjsr(t) 



d m f. 



N 



d(i m 
1 



(12) 



{ki} 



,k l )(\E\ k i)(\E\ k >)...(\E\ k ') 



where the sum over {ki} is taken with the constraint 
ki + &2 + • • • + ki = m and E is the energy of a chain 
of length N + 2 with the first and last sites bound. 
The brackets denote a thermal average. The coefficients 
A(k\, &2, . . . , ki) are combinatorial factors independent of 
the system size N. Since the energy at each site is either 
— 1 (bound site) or (unbound site), one has on the real 
t axis (\E\ k ) < N k . Thus 



N 



dl3 n 



< 



P 



J2 A (ki,k 2 , 



This inequality may be written as 



d m f. 



5/3" 



(13) 



(14) 



where g(m) is independent of N and [3. We now consider 
the wth derivative of the average free energy T ', 



d(3 m 



= d-rf 2 E 



jv 



N=0 



d m f N 



(15) 
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It is evident from (|14|) that this series converges abso- 
lutely which allows us to interchange the order of the 
derivative and the sum. In fact, on the real t axis the 
mth derivative of T is bounded by 



gmjr 



d(3 r ' 



< (l-p) 2 g(m)Y /P N N m - 1 



N=l 



r(m) 

\np\ m 



(16) 



where T(m) is the Gamma function. Therefore the 
derivative exists on the real t axis and particularly at 
t = 0. This analysis implies that the average free energy 
is diffcrentiable to all order in (3 at t = 0. 

An insight of the singularity of the free energy may be 
obtained by considering only the poles closest to t = 
in /aKO f° r calculating T. For example for c > 2 these 
poles are located at t = ±ib/N , where b is a constant. 
Approximating the sum on N in (|^) by an integral and 
after rescaling t one finds 



1 r°° 

Tot-- e- x \n{x 2 t 2 + l)dx 
P Jo 



(17) 



This integral is singular at t = as expected. 

So far we analyzed the limit v — > oo, where no unbind- 
ing takes place. Let us now consider finite v with d»1. 
In this case the model is expected to exhibit an unbinding 
transition at Tm — 0(v), in addition to the Griffiths sin- 
gularity at Tq. This may be seen by noting that at low 
temperatures, T > Tg, the two polymers are expected 
to be bound even for finite but large v. The reason is 
that the free energy difference between the state which 
is bound at all strongly binding sites and the completely 
unbound state is 



AF = —pvL + TpLu In s 



(18) 



where u is a constant of 0(1), Ins is the entropy density 
of a random walk and L is the polymer length. Clearly 
the second, entropic, term is negligibly small at T C 
v. Thus the polymers are expected to be bound at low 
temperatures with an unbinding transition taking place 
at T M =0(v). 

The fact that v is finite but large is not expected to 
modify the nature of the Griffiths singularity found for 
v — > oo at T — Tq. At finite v some of the strongly 
binding sites may become unbound even at the low tem- 
perature around Tg- However this effect is small and 
is not expected to modify the exponential behavior of 
the domain size distribution taken in (|J). For T > Tg 
where all homogeneous domains are unbound our model 
is closely related to the Tang-Chate model in which the 
melting transition has been analyzed. Their analysis of 
the melting transition should be applicable to the melting 
transition in this model as well. 



In summary, we have demonstrated that the model 
considered in this study exhibits a Griffith singularity in 
the limit of infinitely strong binding sites v — * oo. It is 
argued that this transition is present even for finite but 
large v, where it is followed at a higher temperature by 
an unbinding transition. 

We thank E. Domany, Y. Frishman and D. Kandel for 
helpful discussions. 
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